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Abstract 

A Borel summation scheme of subtracting the perturbative contribution from the average plaque- 
tte is proposed using the bilocal expansion of Borel transform. It is shown that the remnant of the 
average plaquette, after subtraction of the perturbative contribution, scales as a dim-4 condensate. 
A critical review of the existing procedure of renormalon subtraction is presented. 
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I. INTRODUCTION 



An old problem in lattice gauge theory is extracting the gluon condensate from the average 
plaquette, which in pure SU(3) Yang-Mills theory has the formal expansion 

Ptf) = 0-- ^Tr U D > = % + YQ m{< ^ GG)aA + ° {a6) ' (1) 

n=l 

where (3 denotes the lattice coupling and a the lattice spacing. The difficulty of extract- 
ing the gluon condensate is that the average plaquette is dominated by the perturbative 
contribution and it is necessary to subtract it to an accuracy better than one part in 10 4 . 
The perturbative coefficients c n were computed to 10-loop order using the stochastic per- 
turbation theory but this alone does not achieve the required accuracy Therefore, any 
attempt to extract the gluon condensate using the perturbative expansion must involve ex- 
trapolation of the perturbative coefficients to higher orders and, the perturbative expansion 
being asymptotic, proper handling of them. Since the large order behavior of perturbative 
expansion is determined by the renormalon singularity of the Borel transform, a natural 
extrapolation scheme would be based on the renormalon singularity. A program along this 
line was implemented by Burgio et a!., and the authors obtained a surprising result of power 
correction that scales as a dim-2 condensate 0]. This is in contradiction with the operator 
product expansion (OPE) ([1]) that demands the leading power correction scale as a dim-4 
condensate. 

The claim of the dim-2 condensate was since then reexamined by several authors. In 
obtaining the perturbative contribution, Horsley et al. employed an extrapolation scheme 
based on the power law and truncation of the perturbative series at the minimal element [3[ , 
and Rakow used stochastic perturbation with boosted coupling to accelerate convergence 
[H, and Meurice employed extrapolations based on assumed singularity of the plaquette in 
the complex /3-plane as well as the renormalon singularity, with truncation at the minimal 
element 0. All these studies did not see any evidence of a dim-2 condensate but found the 
plaquette data was consistent with a dim-4 condensate. 

To help settle these conflicting views on the dim-2 condensate we present in this paper a 
critical review of the renormalon-based approach of ^ , and reveal a serious flaw in the pro- 
gram of renormalon subtraction, and show that the plaquette data, when properly handled, 
is consistent with a dim-4 condensate. 

Specifically, we shall show that the continuum scheme employed for renormalon subtrac- 
tion in [2] is not at all a scheme where the perturbative coefficients follow a renormalon 
pattern, and therefore the claimed dim-2 condensate is severely contaminated by pertur- 
bative contribution and cannot be interpreted as a power correction. We then introduce a 
renormalon subtraction scheme based on the bilocal expansion of Borel transform, and show 
that the plaquette data can be fitted well by the sum of a dim-4 condensate and the Borel 
summed perturbative contribution. 



II. RENORMALON SUBTRACTION BY MATCHING LARGE ORDER BEHAV- 
IORS 

In this section we give a critical review on the renormalon subtraction procedure of 
0. The perturbative coefficients c n of the average plaquette at large orders are expected 
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to exhibit the large order behavior of the infrared renormalon associated with the gluon 
condensate, but the computed coefficients using stochastic perturbation theory turn out to 
grow much more rapidly than a renormalon behavior. This implies that the coefficients are 
not yet in the asymptotic regime, which is expected to set in around at order n = (3z (z 
given below in ([2])), which gives n ~ 30 for /3 ~ 6, far higher than the computed levels. It 
therefore appears all but impossible to extract the gluon condensate directly from using the 
stochastic perturbation theory, since the perturbative contribution must be subtracted, at 
least, to orders in the asymptotic regime. 

In Ref. [2[ this problem was approached by introducing a continuum scheme in which the 
renormalon contribution is subtracted by matching the large order behavior in the continuum 
scheme to the computed coefficients in the lattice scheme. Specifically, in order to relate 
c n of the lattice scheme with the renormalon behavior the average plaquette is written, 
essentially, as 

P(fi) = P ren (/3 C ) + 6P(p e ) + P NP (/3) , (2) 

where 

with P c denoting the coupling in the continuum scheme defined by 

& = ]9-ri-^ (4) 

and 

16tt 2 204 

Zn = , V = . ( 5 ) 

33 121 w 

In Eq. (j2J) the plaquette is divided into perturbative contributions, comprised of the renor- 
malon contribution P ren and the rest of the perturbative contribution 5P, and nonpertur- 
bative power correction Pnp- In this splitting, the asymptotically divergent behavior of the 
perturbative contribution is contained in p ren ; and SP denotes the rest that can be expressed 
as a convergent series. (Here, the renormalons other than that associated with the gluon 
condensate and the subleading singularities at b = z are ignored, which, if necessary, can 
be incorporated in p rcn .) 



We now define P^ p with 



N 



p$P(P) = p(p) - p ren (/3 c ) - $>„ - c™)r n (6) 



n=l 



where C^ en denotes the perturbative coefficients of P ren in power expansion in 1//3. Note 
that Pj^} is free of perturbative coefficients to order N. The constants r 1; r 2 that define the 
continuum scheme and the normalization constant M are determined so that C^ en converges 
to c n at large orders. In the continuum scheme with 

ri = 3.1, r 2 = 2.0 (7) 
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TABLE I: The perturbative coefficients of the average plaquette in the continuum scheme. 

and an appropriate value for Af it was observed that C^ en converge to c n at the orders 
computed in stochastic perturbation theory. The last term in OH]) being a converging series 

P^P will be well-defined at iV — > oo, and this is precisely the quantity that was assumed 
to represent the power correction, and it WclS -f]\fp t licit was shown to scale clS cl dim-2 
condensate. 

The essence of this procedure is that the isolation of the renormalon contribution is 
obtained by matching the large order behaviors in the lattice and continuum schemes, in 
which the matching does not involve the low order coefficients. Although the renormalon- 
caused large order behaviors of any two schemes can be matched, independently of the low 
order coefficients, it must be noted that the matching would work only when the known 
coefficients in both schemes exhibit renormalon behavior. Since, however, the computed 
coefficients in the lattice scheme are far from being in the asymptotic regime and do not 
follow the renormalon pattern the matching cannot be performed reliably; Therefore, the 
conclusion of a dimension-2 condensate based on it should be reexamined. 

That the above matching has a serious flaw can be easily shown by mapping the pertur- 
bative coefficients in the lattice scheme to the continuum scheme ([7]) . If the latter is indeed a 
good scheme for renormalon subtraction the mapped coefficients should exhibit a renormalon 
behavior. However, as can be seen in Table [H which is obtained by mapping the central 
values of c n from the stochastic perturbation theory, the coefficients are alternating in sign 
and far from being of a renormalon behavior. This shows that when mapping the perturba- 
tive coefficients between the lattice scheme and ([7]) the relatively high order coefficients (say, 
7-10 loop orders) are still very sensitive on the low order coefficients. Therefore, the above 
large order matching cannot be performed reliably with the computed coefficients, and (J7J) 
cannot be the right scheme where one can isolate and subtract the renormalon contribution. 

Checking the internal consistency of the subtraction scheme based on the matching of 
large order behavior also shows the underlying problem. The nonperturbative term in (T5]) 
can be written using ([6]) as 

p NP (/3) = pw^) _ {SP{f3c) _ Y^( Cn - on . (8) 

n=l 

For P^p to represent the power correction it is clear that 



A' 



n=l 



(9) 



must be satisfied. Since 5P((3 C ) is by definition a convergent quantity it can be written in a 
series expansion 

oo 

8P%) = J2 D ^c n , (10) 

n=l 
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where D n can be computed up to the order c n are known, and ([9]) can be written approxi- 
mately as 

Now, in the scheme of (JTj), and at iV = 8 and /3 = 6.0,6.2 and 6.4, for example, the ratios 
are 69, 59 and 42, respectively: a severe violation of the consistency condition. This again 
confirms that (J7J) cannot be a scheme suited for renormalon subtraction. 



III. RENORMALON SUBTRACTION BY BOREL SUMMATION 

It is now clear that one cannot subtract the perturbative contribution in the plaquette 
by mapping the renormalon-based coefficients in a continuum scheme to the lattice scheme, 
and then matching them with the computed high order coefficients. On the other hand, 
the lesson of our review suggests that one must map the known coefficients in the lattice 
scheme to a continuum one and look for a scheme where the mapped coefficients follow a 
renormalon behavior. 

Once such a scheme is found one can perform Borel summation to subtract perturbative 
contribution to isolate the power correction. Borel summation is especially suited for this 
purpose, since it allows a precise definition of the power corrections in OPE The 
nature of the renormalon singularity, hence of the large order behavior of perturbation, 
was obtained through the cancellation of the ambiguities in Borel summation and power 
corrections jsj . An extensive review of renormalons can be found in [To| . 

In this paper we shall assume that such a scheme exists and perform Borel summation 



using the scheme of bilocal expansion of Borel transform [llj. To Borel-sum the divergent 
perturbation to a sufficient accuracy for the extraction of power correction, one must have 
an accurate description of the Borel transform in the domain that contains the origin as well 
as the first renormalon singularity in Borel plane. Bilocal expansion is a scheme of recon- 
structing the Borel transform in this domain, utilizing the known perturbative coefficients 
and properties of the first renormalon singularity. After Borel-summing the perturbative 
contribution the sum of the Borel summation and a dim-4 power correction can be fitted to 
the plaquette data. A good fit would suggest then the power correction be of dim-4 type. 

The Borel summation using the first N-loop perturbations of the plaquette in bilocal 
expansion in a continuum scheme is given in the form: 



p(N) 
-'BR 



(/?)= fV/w E V + 7T 

JO Z^n nl (1 



db, (12) 



where the integration over the renormalon singularity is performed with principal value 
prescription. The essential idea of the bilocal expansion is to interpolate the two perturbative 
expansions about the origin and about the renormalon singularity to rebuild the Borel 
transform. By incorporating the renormalon singularity explicitly in the expansion it can 
extend the applicability of the ordinary weak coupling expansion to beyond the renormalon 
singularity, and this scheme was shown to work well with static inter-quark potential or 



heavy quark pole mass [111, [12] . Here, denotes the normalization constant of the large 



order behavior and the coefficients h n are determined so that the Borel transform in (fT2 
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reproduce the perturbative coefficients in the continuum scheme when expanded at 6 = 0; 
Thus h n depends on the continuum perturbative coefficients as well as N '. By definition, 
•^BR when expanded in 1//3, reproduces the perturbative coefficients of the average 
plaquette to N-loop order that were employed in building the Borel transform. For details 



of the bilocal expansion of Borel transform we refer the reader to [12 . 
The power correction can then be defined by 

P^(P)=P((3)-P££(P), (13) 

which, by definition, has vanishing perturbative expansions to order N. 

Using the perturbation to 10-loop order of the plaquette we compute P^r{P) m f ne 
continuum scheme parameterized by Eq. (j3J). Although Af can be computed perturbatively, 
using the perturbations of the average plaquette, it is still difficult to obtain a reliable result 
using the known coefficients, so here it will be treated as a fitting parameter. Thus in our 
scheme, as in the fitting parameters are Af and p,r2 of Eq. (]3J). 

Using the plaquette data for 6.0 < (5 < 6.8 from Il3j] and the relation between the lattice 
spacing a and [3 from static quark force simulation [14J 

log(o/r ) = -1.6804 - 1.7331(/3 - 6) + 0.7849(/3 - 6) 2 - 0.4428(/3 - 6) 3 (14) 

the fit gives Af = 165 and 

n = 1.611, r 2 = 0.246, (15) 

which values are substantially different from those in (J7J). The result of the fit is shown in 
Fig. [TJ which shows that the power correction is consistent with a dim-4 condensate. The 
agreement improves as (3 increases, albeit with larger uncertainties; The deviation at low (3 
(j3 < 6) may be attributed to a dim-6 condensate, which may be seen, though not presented 
here, by that adding a dim-6 power correction in the fit improves the agreement in the whole 
range of the plot. The error bars are from the uncertainty in the simulated perturbative 
coefficients of the plaquette. The uncertainty in the normalization constant does not appear 
to be large: for example, a variation of 20% in Af causes less than a quarter of those by the 
perturbative coefficients. 

From the fit we obtain a dim-4 power correction of Pnp ~ 1-6 (a/ro) 4 . Because of 
the asymptotic nature of the perturbative series the power correction of the plaquette is 
dependent on the subtraction scheme of the perturbative contribution, and thus our result 
may not be directly compared to those from other subtraction schemes. Nevertheless, it is 
still interesting to observe that the result is roughly consistent with 0.4 (a/r ) 4 of [4] and 
0.7 (a/r ) 4 of j5j. Our result turns out to be a little larger than those estimates; This may 
be partly accounted for by the fact that the existing results were from the fit in the low /3 
range of /3 ^ 6, in which range the data are below our fitted curve. 



IV. SUMMARY 

We have reexamined the claim of dim-2 condensate in the average plaquette, and shown 
that the renormalon subtraction procedure of that gave rise to the dim-2 condensate fails 
consistency checks and cannot be reliably implemented with the known results of stochastic 
perturbation theory. We then introduced a renormalon subtraction scheme based on the 
bilocal expansion of Borel transform and found that the plaquette data is consistent with a 
dim-4 condensate. 
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FIG. 1: logPNP vs. j3. The solid line is for 41og(a/ro) + 0.5. The plot shows the power correction 
should be of dim-4 type. 
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